STOCHASTIC REACTION-DIFFUSION SYSTEMS WITH 
HOLDER CONTINUOUS MULTIPLICATIVE NOISE 



MARKUS C. KUNZE 

Abstract. We prove pathwise uniqueness and strong existence of solutions 
for stochastic reaction-diffusion systems with locally Lipschitz continuous re- 
action term of polynomial growth and Holder continuous multiplicative noise. 
Under additional assumptions on the coefficients, we also prove positivity of 
the solutions. 



1. Introduction 

Reaction-diffusion systems and stochastic perturbations of them play an impor- 
tant role in applications in chemistry, biology and physics |19) . In an abstract 
form, a stochastic reaction-diffusion system can be treated as a stochastic evolution 
equation 

(1.1) dU(t) = [AU{t] + F{U{t))]dt + G(U(t))dW H (t) 

on a Banach space E, which is a space of Revalued functions, defined on a domain 
G. Thus, (jl.ip actually represents a system of r coupled equations. Here, A is a 
diagonal matrix of second order differential operators which describe the diffusion in 
the system and the map F accounts for the reaction in the system and is typically a 
composition operator with components of polynomial growth. The system is driven 
by a cylindrical Wiener process Wh in a suitable Hilbert space H . 

In the case of r = 1, i.e. a single reaction-diffusion equation rather than a system, 
there are many articles concerned with such equations, both in the case of additive 
noise (see [3 [TTl [T2j [8] ) and in the case of (locally) Lipschitz continuous multiplica- 
tive noise (see 01 [16j EH HI])- Stochastic reaction-diffusion systems with locally 
Lipschitz continuous multiplicative noise were considered in [5j. 

In the case where the noise term G is no longer locally Lipschitz continuous the 
techniques from the above references can no longer be used. This is essentially due 
to the fact that for such equations a priori only stochastically weak solutions (or 
martingale solutions) can be constructed. However, stochastic reaction-diffusion 
equations with merely Holder continuous multiplicative noise appear naturally, e.g. 
in scaling limits of interacting particle systems, see e.g. [TH] were for r = 1 noise 
terms with G(u)(x) := \u(x)\i and G(u)(x) = \u(x)(l— u(x))\* l{o<«(s)<i} appear. 

For such stochastic reaction-diffusion systems with Holder continuous multiplica- 
tive noise only few results are available and, to the best of our knowledge, only the 
case r = 1 is treated. In [3J, existence of solutions has been proved under an addi- 
tional boundcdness assumption on the coefficient G. However, a uniqueness result is 
missing in that article, except for the case of locally Lipschitz continuous G. Path- 
wise uniqueness for the stochastic heat equation on R d (i.e. r = 1, A is the Laplace 
operator and F = 0) was proved in |2T] in the case where Wh is replaced with a 
colored noise and G is composition with a 7-Holder continuous function, where the 
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allowed values of 7 depend on the noise. In [20) . the one-dimensional version of the 
stochastic heat equation with white noise was considered, and it was proved that for 
that equation pathwise uniqueness holds for 7 > 3/4. For 7 < 3/4, it was recently 
proved in |17j that solutions are neither pathwise unique nor unique in law. 

In this article, we are concerned with stochastic reaction-diffusion systems on a 
bounded Lipschitz domain 6 C R d of the form 

[Mui(t, x) + fi[x, ui(t, x),... ,u r (t, x))]dt 

t > 0,x G 0,1 = l,...,r 
t>Q,x£d0,l = l,...,r 
x G 77,1 = l,...,r. 

Here, £^ is a uniformly elliptic, second order differential operator in divergence 
form on O and djdvA l is the associated conormal derivative. The functions fi : 
ff x H k — > R arc locally Lipschitz continuous, of polynomial growth and satisfy 
suitable dissipativity assumptions. Typical examples that we have in mind are odd- 
degree polynomials with negative leading coefficients, see Section [6] The functions 
gi.k '■ & x R — > R are locally ^-Holder continuous and of linear growth such that 
the Holder constants and the coefficients in the linear growth are square summable. 
Finally, the fti_k are independent, one-dimensional Brownian motions. By |17) . 
pathwise uniqueness cannot be expected for equations driven by space-time white 
noise with ^-Holder continuous coefficient. However, our assumption allow us to 
consider noise terms of the form G(U(t))RdW , where G is composition with a locally 
^-Holder continuous function of linear growth , R is a Hilbert- Schmidt operator on 
L 2 (ff) satisfying additional assumptions and Wh is space-time white noise, see 
Example 12.41 

We will make our assumptions precise in Section [2] 

Under these assumptions, we prove pathwise uniqueness (Theorem 13.11) of solu- 
tions to equation (|1.2p on the state space E = C({7) r ' . Our proof follows the ideas of 
Yamada and Watanabe [29] , who proved pathwise uniqueness for finite-dimensional 
SDE with ^-Holder continuous coefficients. The main difficulty to extend these 
results to the infinite-dimensional setting is of course to handle the differential op- 
erators s^i. The strategy from [20l [21] to prove pathwise uniqueness, which is also 
an adaption of the Yamada- Watanabe ideas, cannot be used in our situation. In- 
deed, there the authors convolute solutions u of the stochastic heat equation with 
a mollifier ip n . In their variational framework, this yields the term u * Aip n in 
the equation for the resulting process. It is then used that, as a consequence of 
its translation invariance, the Laplacian commutes with convolutions, i..e. we have 
u * Aip n = A (it * ip n ). This is no longer true for differential operators with non- 
constant coefficients, which we consider here. We overcome this difficulty by using 
the concept of an (analytically) weak solution, see Definition 12.71 which allows us 
to perform pointwise estimates. 

We would like to point out that in the Yamada- Watanabe result, it is essential 
that the multiplicative noise is diagonal. Therefore, in (|1.2|) we have also considered 
"diagonal noise" by letting the noise term in the Z-th equation only depend only on 
Ui, rather than the whole vector u. Thus, the equations in (|1.2j) are only coupled 
via the reaction terms fi 

Under additional assumptions on the functions fi and g^k, similar arguments as 
in the proof of Theorem l3. II can be used to show that the solutions to (|1.2p preserve 
positivity, i.e. if £1, . . . , > a.s. then we also have Ui{t) > almost surely for all 
t > and I = 1, . . . , r, see Theorem 14.31 



(1.2) 



dui(t,x) 
. uz(0,x) 
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Using pathwise uniqueness, we can adopt the strategy from [51 [T3] to prove 
existence of solutions to (jl.2[) in Theorem 15.51 In contrast to the existence result 
from [3], we can drop the uniform boundedness assumption on G and allow G of 
linear growth. 

This article is organized as follows. In Section^ we fix our assumptions on equa- 
tion (|1.2[> and rewrite it in the abstract form (jl.ll) . We also prove some preliminary 
results and recall some stochastic concepts that we will use. Section [3] contains the 
prove of pathwise uniqueness and Section 2] our result about positivity of solutions. 
Existence of solutions will be proved in Section [5j In the concluding Section [6l 
we apply our results to a stochastic reaction-diffusion system of Fitzhugh-Nagumo 
type. 

2. Preliminaries 

In this section, we fix our assumptions on equation (|1.2j) . rewrite the equation in 
the abstract form and recall different notions of existence and uniqueness for 
equation (jl.ll) that will be used in what follows. 

Throughout, all vector spaces are real. If H is a Hilbert space, we write (•, -)h for 
the inner product in H. When P n (u) and Q n {u) are certain quantities depending 
on an index n and a function u, we write P n (u) < Q n (u) to indicate that there 
exists a constant c, independent of n and u such that P n (u) < cQ n (u) for all n and 
u. We write P n {u) ~ Q n (u) if both P n (u) < Q n {u) and Q n {u) < P n {u). 

2.1. The differential operators. We assume that the differential operators s^i 
are given by 

d 

where we make the following assumptions: 

(A) The domain G C M, d is bounded and has Lipschitz boundary. For I = 
1, • • ■ , r, the matrix valued functions ai := (aL) : 6 — > M, dxd are symmetric 
and have measurable entries. Moreover, for certain 77, M > we have 

d 

ii\Z\ 2 < E a ^Mo < M\t\ 2 

i,i=i 

for almost all x G 6 and all £ eR d . Finally, q G L°°(<^). 

To construct realizations of the differential operators which are generators of 
strongly continuous semigroups, we follow a variational approach. We consider on 
L 2 (ff) the symmetric form 




endowed with domain D(a ) := H 1 ^). The associated operator A;^ is given by 

D(Ai^) := {u e H\0) : 3w G L 2 (^) s.t. (w,v) L 2 = a l [u,v] for allv G iJ 2 (^)} 

and Ai_2U — —w. Note that the boundary condition du/du^ := a;Vu • // = 
on 9^ is incorporated in the domain of Ai, 2. As is well-known, the operator A^2 
generates a strongly continuous and analytic semigroup S^2, see [26] . Changing the 
function q by a constant if necessary, we may and shall assume that is uniformly 
exponentially stable, hence A\^ is invertible. Note that we can compensate the 
change in q by appropriately changing /; in equation (|1.2p . thus this assumption 
means no loss of generality for the stochastic reaction-diffusion equation (| 1 . 2[> . 

Since the form a is sub-Markovian, see [26j Chapter 4], the semigroup Si^ re- 
stricts to a strongly continuous and analytic semigroup Si tP on L p {0) for every 
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p E [2,oo). The generator Ai p of the restricted semigroup Si„ is exactly the part 
of A l>2 in L p (0). 

Finally, D(A L2 ) C C(ff) and the semigroup S1.2 restricts to a strongly continuous 
semigroup S^c on C{0), whose generator Ai t c is the part of A^2 in C(ff). In the 
case where c/ = 0, this follows from the results of [10], the case of general ci £ L°°(ff) 
follows with a perturbation argument. 

To prove existence of solutions to (|1.2|) in spaces of continuous functions, we will 
need some embedding results for the fractional domain spaces D((— Ai :P ) a ). Often, 
such results are obtained from precise knowledge of D(A^ p ) and interpolation. How- 
ever, this strategy requires more regularity of the coefficients ai and the boundary 
of 0. Here, we follow a different approach. 

Lemma 2.1. Assume (A) and let d < 2p. Then for a € (^-,1), the fractional 
domain space D((— Ai yP ) a ) is continuously embedded into C(ff). 

Proof. Let us first note that since H x (0) <-> V^{G) for d > 3 (and H x (ff) ^ 
L q ((?) for arbitrary q g (1, 00) for d < 2), it follows from [TJ Section 7.3] that the 
semigroups Si iP are ultracontractive, i.e. for p > 2 the operator Si tP (t) maps L V (G) 
into L°°(^') and we have 

\\SL P {t)\\^ {L P(ff )}L ^(ff)) < cb~& V0 < t < 1 

for a certain constant c = c(l,p) > 0. 

The fractional domain space D((— Ai jP ) a ) is the range of the (bounded) operator 
(-Ai tP )~ a , defined by 

1 r°° 

(2.1) (_^ ip) -«/ := __ y e -i Slpij . )fdtj 

endowed with the norm ||(— Ai fP )f\\a[(-Ai p ) a ) '■— II/IIp- We note that since Si tP is 
analytic, Si tP (t)f 6 D(A^ p ) for all f > and / <E L v {0). As a consequence of [ID] . 
D(j4; )P ) C C{ff), hence the integrand takes values in C{0). 

Thus, to prove that (— Ai tP )~ a f E C(0), it suffices to show that the inte- 
gral in (|2.ip exists as a Bochner integral in C{0). Hence we have to show that 
||i a_1 <Sz,p(t)/||oo is integrable on (0, 00). From ultracontractivity we obtain 

r-^pM/lloo < cr-'-^WfWp vo < t < 1 

which is integrable on (0,1) since a > Since Si !P is uniformly exponentially 
stable, we obtain for a certain constant d > that 

P a_ Sj,p(*)/I|oo < ll'S';,p(l)l|j2f(Lp(^),L~(€?))ll'S , ;,p( i - l)l|j5f(Lp(^),L~(<y))||/||p 
< c'e-^H/llp Vt>l 

which is integrable on (1, 00). It thus follows that the integral in (|2.1I) indeed exists 
as a Bochner integral in C{0). Moreover, we have the estimate 

IIM)- a /IU<||/||p=||(-Ap)/llD((-A,. p )«.). 

This finishes the proof. □ 

Let us now introduce the Banach space E and the operator A that will be used 
in the abstract formulation (|1.1|) of equation (|l.lj) . We set _E := (C(€?)) r , endowed 
with the norm := YI1—1 \\uj\\<x>- The operator A will be the diagonal operator 

diag(^4i i c, • • • j A r _c), with domain D(Ai i c) x • • • x D(A r ,c)- Then A generates the 
strongly continuous semigroup Sc ■— diag(S'i i c, • • * > Sr,c)- 

To simplify notation, we will drop the index C from now on and merely write 
A%, . . . , A r and Si, ■ . . , S r in what follows. 
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2.2. The reaction term. Concerning the functions /,, we make the following as- 
sumptions: 

(F) For I = 1, . . . , r, the function /, : x R fc — !• R is given as 

fl(x,si,...,s k ) := hi(x,s t ) + ki(x,s 1 ,...,s k ) 

where the continuous function ki is assumed to be locally Lipschitz contin- 
uous and of linear growth in the variables s\, . . . , s r , uniformly with respect 
to the first, i.e. there are constants ci,c 2 and L m for m £ IN such that 

r 

\ki(x, si,...,s r )\ < ci + c 2 ^ \sj\ M x £ 0, sj e R 

and 

r 

|fc,(x, Si, . . . , S r ) ~ h{x,t\, ■ ■ ■ ,U)\ < L m 2J \Sj ~tj\ Vx <E 0, Sj,tj £ [-TO, m] . 

i=i 

The continuous function /i, : ^ x E E is assumed to be locally Lips- 
chitz continuous in the second variable, uniformly with respect to the first. 
Moreover, we assume that 

(1) For a certain constant a, > and an integer TV,, we have 
-a,(l + \s\ 2N ' +1 l {s > 0} ) < hi(x,s) < oi(l + | S | 2Ari+1 l {s < 0} ) 

for all x £ and s £ R. 

(2) For certain constants ai,;, a 2j , £ R and & 2) , > and the integer TV, 
from (Fl) we have 

oi,l - &us 2JVi+1 < ^(x, s) < a 2 ,, - 6 2 ,,s 2Ar ' +1 • 

Example 2.2. Conditions (Fl) and (F2) are satisfied for functions of the form 

2JV,+1 

hi(x,s) := ^ uij^s 3 

3 = 1 

where iV, is an integer and the coefficients ui j belong to C(0) and the highest order 
coefficients w,,2Ar i+1 satisfy ojl2Ni+i(x) < — e < for some e > and all x £ 0. The 
proof of this fact can be found in Examples 4.2 and 4.5 in |14) . 

We now define the operators Hi : C(0) -» C(<?) by (Hi(u))(x) := hi(x,u{x)) 
and the operators Ki : E —> C(0) by (.Kj(ui, . . . , u r ))(x) := fc,(x, iti(x), . . . , u r (a;)). 
Finally, the map F : E E which appears in (|l.ip will be given by F(u) = 
(Fi(u),--- ,F r (u)), where Fj(tt) = Hi(ui)+K t (u). Obviously, the maps F, F, , if, , if , 
are locally Lipschitz continuous. 

Conditions (Fl) and (F2) imply that the maps H satisfy additional dissipativity 
conditions which play a crucial role in proving existence of solutions of (|l.lj) . These 
dissipativity conditions were first used in [7J. 

Recall, that in a Banach space {B, \\ ■ ||), the subdifferential of the norm at x is 
given by 

:={^eB* : |M| =1 and (x,<p) = ||x||}. 
As a consequence of the Hahn-Banach theorem, d\\x\\ ^ for all x £ B. We have 

Lemma 2.3. Let Hi be defined as above, u,v £ C(0) and ip £ <9||xt||oo- Then we 
have 

(1) (Hi(u + v),<p) < a'(l + |H|oo) 2JVl+1 and 

(2) (H(u + v)- Hi{v),<p) < a"(l + ||Hloo) 2Ar!+1 - b"\\u\\^'+ 1 . 

for suitable constants a',a",b" > depending only on the constants a, from (Fl) 
and au,a 2 ,,,&i,,,6 2 ,, from (F2). 
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Proof. See Section 4.3], cf. also Examples 4.2 and 4.5 in [H] □ 

2.3. The noise term. Let us now turn to the stochastic term in equation (|1.2[) . 
We shall assume: 

(G) For I = 1, • • • , r and k G IN the continuous functions gi_k :#xl-)l have 
the following properties 

(1) There exist a; = (a^ £ K,A = (X)fceiN € ^ 2 such that \gi tk (x,s)\ < 
a[ + (3 l k \s\ for all x G ^ and set 

(2) For every ui (E IN th.6rc exist continuous functions <j\ ^k,m 

: [0, oo) — > 

[0, oo) such that 

(a) We have \gi, k (x, si) - gi )k (x, s 2 )\ < cri, k ,m(\si - s 2 \) for all x G ^ 
and si,S2 G [— m, m]. 

(b) For I = 1, ■ • • , r and m G IN, the function : (0, oo) — > (0, oo), 
defined as p;, m (s) := Y^kLi Ci.fc.mt^) 2 is increasing in s and sat- 
isfies 

/ 7~rds = oo . 

Jo+ piM s ) 

Example 2.4. Let us give an example of functions gi ik which satisfy (G). For sim- 
plicity, we assume that r = 1 and drop the subscript /. 

Let R : L 2 (&) — > L 2 (&) be a Hilbert-Schmidt operator which is diagonalized by 
an orthonormal basis (ek)k£Ki i.e. for some sequence A = (Xk)k&N G we have 

oo 

Ru = ^2\k(u, e fc ) L 2e fc . 

k=\ 

We additionally assume that e k G C(^) and that the sequence || Afcefc||oo is square 
summable. Moreover, let g : R — >• R be of linear growth and locally Holder contin- 
uous of order i, i.e. \g(s)\ < a + b\s\ for certain a, b G R and there exist constants 

c m > such that |g(si) — <?(s2)| < Cm|si — «2 1 ^ for all si,s 2 G [— m, m]. 

Then (G) is satisfied if we set gk(x, s) := g{s)Xkek{x). Indeed, (1) holds true with 
a k = allAfeefeHoo and /3 k = ^Afce/s^oo. For (2), we choose a k , m (s) = c m \\X k e k \\ooS^ ■ 
Then clearly (a) is satisfied. In (b), we find p m (s) = sJ^kLi 1 1 A/c e/c 1 1 ^3 , hence p m x is 
not integrable near 0. 

In this situation, the noise term in equation (|1.2j) can actually be rewritten as 
g(u(t,x))RdW(t,x), where W is a space-time white noise, i.e. an L 2 (i£?)-cylindrical 
Brownian motion. 

Remark 2.5. If the orthonormal basis (efc) in example 12.41 is uniformly bounded in 
C{&), i.e. sup fc ||efe||oo < 00, then our assumptions reduce to (Afc) G i 2 . This is for 
example the case if we consider the standard orthonormal basis of L 2 (0, 2tt). If we 
want to consider the orthonormal basis (e k ) which diagonalizes the operator A 2 on 
L 2 (i^), then this basis consists of functions in C(ff), as D(A S ) C C(i?), however the 
functions are typically not uniformly bounded in C{ff). Typically, one can obtain 
a bound IjefcHoo < k s for a certain s > from ultracontractivity. In this case, one 
has to color the noise more by requiring that (k s Xk) G £ 2 . 

Let us return to our general setting and explain how the stochastic term is 
modeled in the abstract equation ([TTT]) . We define G t : C{G) ^{i 2 ,C(ff)) by 

OO 

[Gi{u)h] (x) :=J2h k [Gi, k (u)](x), 

k=l 

where G itk : C{ff) C(W) is defined by [Gi, k (u)](x) := gi <k (x,u(x)). Note that 
the above series converges in C(^), since \\h k gi lk (-, u (-))lloo < h k (at tk + A,fe|| u ||oo) 
is summable. 
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For the purpose of stochastic integration, we will view G; as a function taking 
values in ^((. 2 ,L P {0)) by embedding C(G) into L p (d). It turns out that Gi even 
takes values in -y(£ 2 , L p (ff)), the space of j-radonifying operators from I 2 to L v (6\ 

Let us briefly recall the definition of 7-radonifying operators. Given a Hilbert 
space and a Banach space B, every finite rank operators T : J4? — > B can be 
represented in the form 

A 

T = h j ® Xj 
i'=l 

for some integer iV > 1, where (hj)jL-y is an orthonormal system in and 
afi, . . . , Xjv 6B. Here, /i (g> x is the rank one operator mapping 5 G J^ 5 to (5, h)^>x. 
With T represented in this form, we define 

A „. 1 



\\T\\y(je,B) ■= { E \\}_^lj x . 

j=i 

where (7j)|Li * s a sequence of independent standard normal random variables. 
This norm is independent of the representation of T in the above form. The space 
"f(Ji?,B) is the completion of the finite rank operators with respect to the norm 
II ' W-y(^f.B)- The identity operator on Jrff <g> B extends to a continuous embed- 
ding 7(J0°, B) J£{J%',B). We may thus view j(Jif,B) as a linear subspace of 
_Sf (J^, B). The operators belonging to "fiJ^f, B) are called the 7-radonifying oper- 
ators from Ji? to B. For more information about 7-radonifying operators, we refer 
the reader to [25] . 

To see that the map Gi from above takes values in 7(f 2 , L p (&)), let (e^) denote 
the canonical orthonormal basis of £ 2 . Then the series 



£|G*(n)e fe | 2 = X>.fcM-))| 2 
fc=i fe=i 

converges in W{G). Thus, by [25] Lemma 2.1], G{u) G 1 {l 2 ,W{e)) and 

OO l 

l|G(«)[[ 7(<a , w(<T) ) ^ |(Ek*(-.«(-))| a ) , || " 1^1(11"'^ + INMMloo), 



<?2 



proving that Gi is of linear growth. Let us next show that Gi : C{0) — > j(£ , L p (ff)) 
is continuous. To that end, let u n — > u in C{&). Employing [25j Lemma 2.1] a 
second time, it follows that 



\\G(u n ) - G(u)\\ y (p,Lr{0)) ^ lflU0> "»(*)) -5/,fc(-,w(-))| : 

fe=i 

which converges to as n — > 00 by dominated convergence. 
Let us summarize the properties of Gi 

Lemma 2.6. Assume (G). Then the maps G t : C{0) -» ^{i 2 ,LP{&)) are well- 
defined, of linear growth and continuous. 

We now proceed to model the stochastic term in equation (|l.lj) . The driving 
process Wh is an ff-cylindrical Wiener process for a suitable Hilbert space H, 
defined on a filtered probability space (0,E,F,P), i.e. a bounded linear operator 
from L 2 ([0, 00); H) to L 2 (Vl) with the following properties: 

(1) for / G L 2 ([0, 00); H), the random variable Wn(f) is centered Gaussian. 

(2) for t > and / G L 2 ([0, 00); i?) with support in [0,t], W H (f) is 
measurable. 

(3) for t > and / G £ 2 ([0, 00); H) with support in [t, 00), Wfr(/) is indepen- 
dent of 

(4) for h,f 2 G i 2 ([0,oo);i?) we have E(W H (fi)W H (f 2 )) = (/1, / 2 )^([o l0 o) ; ff). 
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It is easy to see that for h € H, the process (Wh {t)h)t>o, defined by Wn{t)h := 
Wij(l(o,t] 8> h) is a real valued F-Brownian motion (which is standard if \\h\\jj — 
1). Moreover, two such Brownian motions (Wff(i)/ii)t>o and (WH(t)h 2 )t>o are 
independent if and only hi and h 2 are orthogonal in H . We refer to [23 for a 
further discussion. 

To model the reaction-diffusion system (|1.2|) in abstract form, we choose i? := 
[^ 2 ] . Denoting the canonical orthonormal basis of £ 2 by (e^), we put e;^ := 
(0, • • • ,0, efc, 0, . . . , 0), where the e k is at position I. Then (e^k) is an orthonormal 
basis of H. Let f3^ k , I — l,...,r and fc G IN be a family of independent real- 
valued Brownian motions defined on a common probability space (f2,E,P). Then 
W H : L 2 (R + ,H) -4 L 2 (0), defined by 

r oo „oo 

: = E E / (/(*)> ^MA.fcW 
i=i k=i Jo 

is an if-cylindrical Wiener process, see [23 . We denote by Pi : H —> H the 
projection onto the l-th. component and define G : E — > Jz?(H, E) by 

G(u)/i := (Gi(ui)Pi/i,...,G r (u r )P r /i). 

Later on, we will also write Hi := PiH and define Wh l by Wn^f) '■= WuiPif) for 
/ G L 2 ([0, oo); iJ). Then Wjj, is an iJ;-cylindrical Wiener process. 

2.4. Solution concepts. Having rewritten equation (|1.2[) in the abstract form 
(II. ip . we now define what we mean by a solution of equation thus by a 

solution of equation (|1.2[) . Since the map G is not locally Lipschitz continuous, we 
initially work with stochastically weak solutions, i.e. the probability space is part of 
the solution. We first consider a solution concept which is also analytically weak. 

Definition 2.7. A weak solution of equation is a tupel ((f2, E, P), F, Wh, u), 
where (fi,E,P) is a probability space endowed with a filtration F which satisfies 
the usual conditions, Wh is an ii-cylindrical Wiener process with respect to F 
and u = (u(t)) t > is a continuous, F-progressive, .©-valued process such that for all 
x* e D(A*) andt > we have 
(2-2) 

(u(t),x*) = (u{0),x*}+ I (u(s),A*x*) + (F(u(s)),x*)ds+ [ G(u(s))*x* dW H (s) . 
Jo Jo 

P-almost surely. 

If an initial datum £ £ E is specified, we say that ((f2, E, P), F, Wh, U) is a weak 
solution to the initial value problem corresponding to Ijl.ljl , if it is a weak solution 
of (d3J and additionally P([/(0) = £) = 1. 

Note that the stochastic integral in (|2 .2[) is an integral of an H ~ iJ*-valued 
stochastic process with respect to a cylindrical Wiener process. Such an integral is 
defined as follows. If (hf.)kevi 1S an orthonormal Basis of H, put Wk '■= Wn(s)hk 
which is a real valued Brownian motion. We then define 

t oo „f 

G(u(s)Yx*dW H (s):=y2 / (G(u( S )yx*,h k ) H dw k (s). 
k=i J o 

In the special situation of equation (|1.2j) . D(A*) can be identified with D(yl*) x • • • x 
D(A*) which is a subset of ^(ff) r , the r-fold product of the bounded Borel measures 
on G. Using the standard basis (ei t k) of H = [£ 2 ] and writing fak := Wh&u as in 
the previous section, we have for x* = . . . , /x r ) G D(^4*) 

ft r oo „ 

/ G(u(s))*x* dW H (s) = EE _9i,k(x,ui(s,x))dni(x)d(3^ k (s) 
Jo i=ik=i Jo 
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In particular, choosing p,\ :— R(X, A{)*8 X , where S x is Dirac measure in x and fij = 
for j ^ I, equation (|2.2I) reduces to 

[R{X,Ai)ui(t)](x) = ui(0,x)+ [ [A l R(\,A l )u l (s)](x)+[Fi(u( s ))](x)ds 

Jo 

r 

+ Y, [R(X,Ai)Gi(u(s))] (x)dW Hl (s). 
1=1 

This will be used in the following section. 

In order to prove existence of solutions, the concept of weak solutions is not 
suitable. Instead, we will use the concept of a mild solution. To define the concept 
of a mild solution, we have to use a Banach space valued stochastic integral. We 
will use the theory of stochastic integration in UMD Banach space [24. We note 
that our state space E is not a UMD Banach space. However, the fractional domain 
space D((—Ai >p ) a ), being isomorphic to L p {i?), is a UMD Banach space. We may 
thus perform stochastic integration in fractional domain spaces and then use Lemma 
12.11 to get back into our state space. More precisely, we have 

Lemma 2.8. Assume (A) and (G). Moreover, let (fi, E, F, P) be a stochastic basis 
on which an H- cylindrical Wiener process is defined. Then, if X = (X\, . . . , X r ) 
is a continuous, F '-progressive E-valued process, then for every t > the process 
t M> S(t — s)G(X(s)) is stochastically integrable and the stochastic integral 

f S(t-s)G(X(s))dW H (s) 
Jo 

defines an ^ t -measurable, E-valued random variable. 

Proof. It suffices to prove the assertion for r = 1. We will drop the index / for ease 
of notation. Fix uj e fl, t > 0, p > max{2, d/4} and a £ (^, |). We claim that the 
map Ru : L 2 (0,t;H) -> D((~A p ) a ) given by 

(R^f,x*) := / (S(t-8)G(X(8,u)),x*)da 
Jo 

for all x* € D((—Ap) a )* is 7-radonifying. Since D((— A p ) a ), being isomorphic to 
L p (r2), has type 2, it suffices to show that s H> S(t — s)G(X(s,uj)) belongs to 
L 2 (0,t;~/{H, D{{-A p ) a )), see [23 Theorem 11.6]. To that end, we have 

\\S(t - -)G(X(-, w))|| i 2( ,t; T (ff,D((-A p )<*)) 

= f \\S(t - s)G(X( S ,oj))\\* {HtDa _ Ap)a)) ds 

J 

< l\t - s)- 2a \\G(X(s, •))||^.d((-a p )*)) ds 
Jo 

< fit - S )- 2 «(H|H|, 2 + ||/3|H|X(^)||oc) 2 ds < 00, 
Jo 

since a < i. Here, we have used the estimate ||5 , (i)||^f(Lp(^),D((-A p )")) ^5 t~ a and 
the ideal property for 7-radonifying operators 1231 Theorem 6.2] in the second line 
and the linear growth of G from Lemma 12.61 in the third. 

It now follows from Lemma 2.7 and Remark 2.8], that s 1— > S(t — s)G(X(s, lu)) 
represents a strongly measurable map R : O — > 7(L 2 (0, t; H); D((— A p ) a )). As 
s i-> S(t — s)G(X(s,ui)) is F-progressive, it follows from Theorem 5.9] that the 
process is stochastically integrable and the stochastic integral defines a D((— A p ) a )- 
valued random variable. By the embedding from Lemma |2.1[ we are done. □ 

We may thus define 
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Definition 2.9. Assume (A), (F) and (G). A mild solution of equation (jl.lj) is 
a tupel ((f2, E, P), F, Wh, u), where (Q, E, P) is a probability space endowed with 
a filtration F which satisfies the usual conditions, Wh is an f/-cylindrical Wiener 
process with respect to F and u — (u(t)) t > is a continuous, F-progressive, i5-valued 
process such that for all t > we have 



F-almost surely. 

If an initial datum (, E E is specified, we say that ((O, E, P), F, Wh, U) is a mild 
solution to the initial value problem corresponding to (ll.ip , if it is a mild solution 
of CLU) and additionally P(£/(0) = £) = 1. 

By the results of 13, Section 6], the weak and the mild solutions of equation 
(II. ip coincide. We will thus briefly speak of solutions of equation (II. ip . rather than 
weak (or mild) solutions. 

2.5. Pathwise uniqueness and strong existence. Typically, when working with 
stochastically weak solutions, the appropriate uniqueness concept is that of unique- 
ness in law. However, we will use the following uniqueness concept: 

Definition 2.10. We say that pathwise uniqueness holds for (11.11) . if whenever 
((f2, E, P), F, Wh , Uj) for j = 1,2 are weak solutions of equation (|l.ip . defined 
on the same probability space and with respect to the same 77-cylindrical Wiener 
process, satisfying iti(O) = 1/2(0) a.s. we have P(ui = U2) = 1. 

For finite dimensional stochastic differential equations, Yamada and Watanabe 
[29] proved that pathwise uniqueness implies uniqueness in law. Moreover, they 
proved that if solutions exist, then they exist strongly, i.e. given a stochastic basis 
(O, E,F,P) on which an iJ-cylindrical Wiener process is defined, we can find a 
continuous, F-progressive, -E-valued process u defined on that stochastic basis, such 
that ((fi, E, P), F, Wh, u) is a weak solution. 

These results also generalize to stochastic equations on Banach spaces, see |T3l 
Section 5]. 

We will make extensive use of strong existence of solutions in Section [5] to prove 
existence of solutions to stochastic reaction-diffusion systems also for unbounded 
reaction terms. 

We should also note that given pathwise uniqueness and existence of solutions 
for deterministic initial values £, if follows automatically that solutions exist for 
random initial data £ : il —> E which are J^o-measurable. See [13] for a proof in the 
infinite dimensional case. 



In this section we prove 

Theorem 3.1. Assume (A), (F) and (G). Then pathwise uniqueness holds for 
equation (|1.2p 

Proof. We may (and shall) assume that pi, m (t) > t for I = 1, . . . , r and m € IN, 
otherwise replacing pi, m {t) with pi, m {t) + t. 

Let us fix I G {1, • • • ,r} and m S IN. Similar as in the proof of the classical 
Yamada- Watanabe theorem [29] . we chose a decreasing sequence a n I such that 
ao = 1 and 




3. Pathwise uniqueness 
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This is possible by assumption (G2b). Note that the sequence (a n ), as well as the 
functions ip n and tp n introduced next, depend on I and m. Next we pick functions 
ipn € C£°(It) such that suppi/Vi C (a n ,a n -i) and 

< Mt) < ~~77\ < and / i> n {t) dt = 1. 

npi t m(t) nt J an 

We define 

/•|t| 

¥>n(t) := / / ip n (r) dr ds . 



Then <p n £ C°°(R) with 

rl*l 

= s § n (*) / ^n(s)ds and = ^„(|t|). 

Jo 

We note that 

,1*1 

|*| -On-l = / l(a n _i,oo)(s)ds < < 

Jo 



which implies that i/3 n (t) — > \t\, uniformly on It. Moreover, f' n (t)t = \t\ jj*' ip n (s) ds 
converges to |i| pointwise. 

After this preparation, we now come to the main part of the proof. Let u\ :— 
(u\,..., u].) and U2 — (uf , . . . ,v%) be two weak solutions of (|1.2p . defined on the 
same probability space and with respect to the same sequence of Brownian motions 
(fil t k)- Moreover, we assume that iti(O) = 112(0) almost surely. 

For m £ IN we define the stopping time r m by 

r m := inf {t > : |M*)|| B V \\u 2 (t)\\ E > m}, 

where we set inf = 00. 

For A > 0, I £ {1, . . . , r} and x £ & the vector x* — (/ii, . . . , fi r ), defined by 
fi* := XR(X,Ai)*5 x and = for j ^ I, is an element of D(A*). Thus, since u\ 
and M2 are weak solutions with Mi(0) = 1*2(0) almost surely, we have, almost surely, 

[XR(X,Al)(ul(t A r m ) — uf(t A r m )] (x) 

[AiXR(X,A l )(uj(s)-uf{s))] (x)ds 

+ / [XR(X,Ai)(F l (u 1 ( s ))-F l (u 2 ( s )))](x)ds 
Jo 

+ E / l[o.r m] (s)[Ai?(A,A i )(G i , fc ( Mi 1 ( S ))-G / , fc (u2( s )))](x)dA, fc ( S ). 

fcl- 70 

To simplify notation, we introduce some abbreviations. We write 
A x ui(t) := XR(X,Ai)[uj(tAT m )-u?(tAT m )] 
AxFtit) := XR{X,Ai)[{ui{tAT m ))-Fi{u2{tAr m ))] 
&xGi,k{t) ■= ^(A,i4i)[Gi )fc («f(tAT m ))-Gi lfc («?(tAT m ))]. 
From Ito's formula it follows that, almost surely, 

<p n (A x ui{t){x)) = / ^(A A m(s)(a;))[AjA A u z (a)(a;) + A x F l {s)(x)] ds 
Jo 

+ -/ ^(AaU{(s)(x)) ^ [A\Gk,i(s)(x)] ds + a martingale. 

fc=i 
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As this is true for every x 6 we may integrate over and take expectations. 
This yields 

E (p n (A xUl {t){x))dx = E / / (p' n (A x u l (s)(x))A l A x ui(s)(x)dxds 



+ E / ^(A A ^(s)(x))A A F ; (s)(z)cfeds 

(3.1) 



^(A A «,(s)(a:)) ^ [A A G;, fe (s)(z)] dx ds 



fc=i 

=: / 1 (n,A) + /2(n,A) + 7 3 (n,A). 
We proceed in several steps. 

Step 1: We estimate A) and let A — >• oo. 

Noting that ^(A A u;(s)) S H 1 ^) by the chain rule, we find that 



7i(n,A) = E (A i A A u ; ( s ),^ l (A A u z (s))) i2(€?) 



ds 



thr„ 



E a l [A x ui(a),(f/ n {A x ui(a))]d8. 



Next observe that 

a 1 [A x ui(s),<p' n (A x ui(s))] 

P d 

= 12 a io ■ (AAa«i(»)) • ^(A A «j(a)) • (£> i A AU ;( s )) dx 

d 

> / r)^(A x ui(s))Y\ \D t A x ui(s)\ 2 dx > , 

since ip'^ > 0. It thus follows that Ii(n, A) < 0. 
We now abbreviate 

Aui(t) := u}(t A T m ) - uf{t A r m ), AFi(t) := Fi{ui(t A r m )) - F ; (u 2 (i A r m )) 
and 

AGi lfc (t) := Gi, k {ul(t A r m )) - Gj, fc (u?(t A r m )). 

For w e C(£?) we have that Ai?(A, converges to u> in C(ff) as A — >■ oo. Thus, 
it follows that A x ui(t) -> Au/(i), A x Fi(t) ->• Af;(t) and A x Gi t k(t) -> AG/,fe(*) in 
C(i^) as A — ► oo for every i > 0. Inserting the estimate for Ii(n, A) into (|3.1|) and 
letting A — > oo we obtain, using the continuity of ip n , cp' n and 

E ip n (Aui(t)(x))dx <E / / ^(Aui(s)(a;))AFi(s)(a;)rfa;rfs 
JV Jo JV 



-E/ / </?" (Aui (s)(x))^2 [^ G l,k (s) (x)] 2 dx ds 
2 ^ Jff fe=i 



2 

=: Ji(n) + J 2 (n) 
5iep ,2: We estimate Ji(n) and J2(n) and let n — > oo. 

Since is locally Lipschitz continuous by (F) and since \ip' n \ < 1, it follows that, 
for a constant L rn > 0, we have 

Ji(n)<L m E y2\Auj(s)(x)\dxds < L m E / V |Au j (s)(a;)| dxds . 

Jo ^ ° j=\ J° e j=i 
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As for J-zin), using (G2) and the estimate (p' r [(t) — ip n (t) < — - ^ , we see that 

J 2 (n < o E / / 1 R \ {0 } Am(s) a: ^ fc " 1 V, A , w\ dxds 

K [ tATm [ t, rA " 2p,, m (A( a )(x)) ( fa; t 

= -E / / 1r\ {0 } Au;(s) x T7- — rfs < - 

2 Jo Je npL m {\Aui{s){x)\) n 

Combining these estimates, we find that 



E [ <p n (Aui(t)(x))dx < L m [ [ \Auj(s)(x)\ dxds + - 
Je Jo J&j=i n 

Since <p n (t) t I* I as n — » oo, it follows upon n —} oo that 

(3.2) E [ \Aui(t,x)\dx < L m E [ [ V \Auj{s)(x)\ dxds 

Je Jo J e - =1 



Step 3: We finish the proof. 

As equation (|3.2[) is true for every £ = 1, . . . , r, we find, summing up, that 

E ( y2\Aui(t,x)\dx<rL m [ E [ *^T\Au{s)(x)\dxds 
J e , =1 Jo Je l=1 

Thus, by Gronwall's Lemma, 

E / y~] \Aui(t, x)\dx = E / V \u](t AT m ,x) - uf(t AT m ,x)\dx = 0. 
Je l=1 Je l=1 

Since t A r m —¥ t almost surely, upon m — > oo it follows that 



E f V - uf(t,x)\dx = 

^ i=i 



for all i > 0. As solutions are continuous in x, it follows that u\(t) = v,2(t) almost 
surely, for every t > 0. Finally, by continuity of the paths, the exceptional set may be 
chosen independently of t, hence u\ — U2 almost surely, i.e. pathwise uniqueness. □ 

Remark 3.2. Note that in the proof of Theorem 13.11 neither the special structure of 
the functions /; in condition (F) , nor the linear growth assumption on the functions 
gi t k in condition (G2) was used. Thus pathwise uniqueness holds already if the 
functions /; are locally Lipschitz continuous and the functions gi t k satisfy (G2). 



4. POSITIVITY 

Under additional assumptions on the nonlinearities /; and g^i, the techniques 
from Section [3] can also be used to prove that the solutions of equation (II. 2[) for 
positive initial data are almost surely positive. For deterministic reaction-diffusion 
equations, the solutions preserve positivity if and only if the reaction term is quasi 
positive [28] . 

Definition 4.1. A map $ = (<£>i, . . . , $ r ) : 6 x W — > W is called quasi positive if 
for every I = 1, . . . , r and s — (si, . . . , s r ) with si — and Sj > for all j =/= I wc 
have $i(x, s) > for all x e ~& . 

In what follows, we write s + := s V for the positive part of a real number s and 
s~ := (— s) V for the negative part of s. 



14 



MARKUS C. KUNZE 



Lemma 4.2. Let $ : ffxR k -» R fe be quasi positive and locally Lipschitz continuous 
in the sense that for ngl there exists a constant L n such that 

r 

\$l(x,si, ...,s r )-$i(x,ti,.. .,t n )\ < L n ^2 \ s j — *j I ^ x S &,Sj,tj G [-n,n] . 

3 = 1 

Then for every I G {1, . . . , r}, x G <*? and s\, . . . , s r G [— n, n] with si < we have 

r 

(4.1) -<$>t(x,si,...,s r ) < L n ^2sJ. 

3=1 

Proof. Without loss of generality, we assume that Z = 1. We have 
-$i(a;, si, . . . , Sr) = $i(x,0, s 2 , . . . , s r ) - $i(x, si, . . . , s„) - 0, s 2 , ■ ■ ■ , s r ) 

< L„|0-Si| -$> 1 (x,0,s 2 ,...,s r ) 

= L„S^ - $i(x,0, S 2 , . . . , S r ) . 

If <I>i(x, 0, s 2 , . . . , s r ) > 0, this already yields (I4.ip . Otherwise, since $ is quasi 
positive, there exists an index 2 < jo < r with Sj < 0. We assume without loss of 
generality that jo = 2. The same estimate as above shows that 

0,S 2 , . . . ,Sr) < L n Sz - $i(x,0, 0,S 3 , . . . ,S r ) . 

These arguments are now iterated. Doing this, we have proved (|4.1I) at some point 
or, after r iterations, we have proved that 

r 

-$i(x,si, . . . ,s r ) < L n ^2sJ - $i(x,0, ... ,0) 

3=1 

which also yields (|4.1j) . as $i(x, 0, . . . , 0) > by quasi positivity. □ 

Theorem 4.3. Assume (A), (F) and (G). Additionally, assume that f : & x W — > 
MT is quasi positive and that gi,k{x, 0) = for all x G G , I = 1, . . . , r and tel. 

If u = (ui, . . . , M r ) is a solution of (|1.2[) im'i/i ui(0), . . . , u r (0) > almost surely, 
then almost surely also ui(t, x) > /or all t > and x G . 



Proof. We pick the sequence a„ and the function ip n as in the proof of Theorem 13.11 
and define ^ 

tp n (t) := l(o,oo) / / ip n (r)dTds. 
Jo Jo 

Then <^„ G C°°(R) with t £ + - Moreover, 

= ^(O.oo) / ipn(s)ds and <^"(t) = l (o , oo) -0 n (i) . 
Jo 

In particular, </?^(i) t l(o,oo) f° r an i G R. 

Given a solution u of (jl.2p with ui(0), . . . , « r (0) > almost surely, we repeat the 
computations and estimates in the proof of Theorem 13 . 1 1 with u\ = and m 2 =: u. 
Doing so, we obtain at the end of Step 1 

tp n (-ui(t A T m ,x)) 

Iff 

rfAT m p 

< E / / ip' n (— ui(s AT m ,x)) ■ (-fi(x,u(s Ar m ,x))dxds 
Jo Jff 

+ -E / / <p"(-ui(s Ar m ,3;))y [-9i,k(x, ui(s Ar m ,ac))] dxds 
2 ^ fc=i 
=: Ji{n) + J 2 {n). 

Note that we do not obtain an extra term for the initial datum, as ip n (—vi(x)) = 
almost surely. We now estimate the terms J\ (n) and J 2 (n) . 
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Since <p' n (t) — for t < 0, the integrand in J\{n) vanishes, unless u;(sAr„, x) < 0. 
In that case, since / is quasipositive and locally Lipschitz continuous, it follows from 
Lemma 14.21 that 

Ji{n)<E / ip' n (-ui(s A r m , x))L m Uj(s A r m , x)~~ dx ds . 

Jo Je j=1 

Concerning J2(n), since gi,k{x,ti) = 0, we obtain 

= -E / / <p"(~ui(s Ar m ,x)) V] [gi,k{x, 0) - gi.k{x, u t (s Ar m ,x))] cteds 

. 1„ f tArm /" 2pi(—ui(s A T m , x)) t 

< -E / / l(o l0 o)(-«i( s A 7>n,x)) ? 1— --dxds<- 

2 io Jff npi(-ui{s AT m ,x)) n 

With these estimates and n — > oo, it follows that 

E / ui(tAT m ,x) dx < L m I E / N Uj(s Ar m ,x) dxds 



e 



f E f £u,(* 
Jo Je =1 



since ip n (-ui(t A r m ),x)) t ui(t A r m , x)~ . 

Now the same arguments as in the proof of Theorem 13.11 show that Y%=i u i = 
almost surely, hence, almost surely, ui > for all I = 1, . . . , r. □ 

5. Existence of solutions 



We now prove existence of solutions to equation (|1.2|) . In the proof, we will 
use the pathwise uniqueness proved in Section |3l Indeed, by pathwise uniqueness, 
solutions exist strongly, i.e. on a given probability space and with respect to a given 
cylindrical Wiener process. This allows us to adopt the strategy of [2] to our 
situation and prove existence of solutions via approximation of the coefficients. 

We will use the following estimates for deterministic and stochastic convolutions 
which are a consequence of the factorization technique [6] . 

Lemma 5.1. Assume (A), (F) and (G). Let (f2,S,F, P) be a stochastic basis on 
which an H -cylindrical Wiener process is defined. Moreover, let u = (u±, . . . , u r ) be 
a continuous, IF '-progressive E -valued process. Then for p > 2 there exists a constant 
C(T) with C{T) — > asT — > 0, independent of the constants c\ and ci from (F) and 
the sequences ai and fii from (Gl) such that for the stochastic convolution process 

SiOiGi{ui) :=t^ f Si(t - s)Gi(ui(s) dW Hl (s) 
Jo 

we have the estimate 

E \\Si o, G l {u l w c{[wm < C(T)(|HI? 2 + BAIW".^,^)) 

and for the deterministic convolution process 

S^Kiiu) :=t4 f Si{t-a)Ki(u(s))da 
Jo 

we have the estimate 

E\\st * < c(t)(c? + 4nn p E ) . 

Proof. We first treat the stochastic convolution. We pick q so large that ^ < \ — ^ 
and a 6 i — |). Then, for some e > and f3 e (a + |, |), we find 

< MSi<>iGi{ui)\\ p c{[0T] . D{{ _ Aig)a) 
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£ ^j[ T E|| a .->(t-«)-^G I («i(*))ll^ 3 (o 1 t^ l ),w W )*' 

< r ei, jr T E(jr*(t- a )- a /'iiG I (t»,)ii» (Hl , i4W) da)*dt 

< T ep / T /"'(f _ s )-2/3|^|2 ( || a; || £2 + ||/3 ; || £2 || U; || c([0iT . c(<y)) ) 2 d S ) 5 dt 



< / E|^(||aj||^ + ||A||/»||«i||c([o,TiOW)) p <ft 

JO 



Here, the first estimate follows from the embedding of Lemma 12. 1[ the second 
is [25j Proposition 4.2], the third uses that L q (0) has type 2, whence we have 
the embedding L 2 (0, t; y(Hi,L q (0)) 7(£ 2 (0, t; Hi), L q (0)). The fourth estimate 
follows from the linear growth of Gi, proved in Lemma \2 .6[ the fifth from the fact 
that sup tg [ T ] f Q (t — s)~ 2 ^ ds < oo and the last one from Young's inequality. This 
proves the first assertion. 

The (easier) proof of the section assertion is similar, using [15, Proposition 4.2.1] 
instead of [25l Proposition 4.2]. □ 

We first prove existence of solutions under additional boundedness assumptions 
on the maps ki and gi^. 

Theorem 5.2. Assume (A), (F) and (G). Additionally assume that 

(1) The functions ki from assumption F are uniformly bounded, i.e. there exists 
a constant C > such that \ki(x,s)\ < C for all (x,s) £ x R k and 
1 = l,...,r. 

(2) For I = 1, . . . , r, the vector fti from assumption (Gl) satisfies Pi = 0. 

Then every £ = . . . ,£ r ) £ E, there exists a unique solution of equation (II. 2ft 
with initial datum £. 

In the proof of Theorem l5.2[ we use the following Lemma which is a reformulation 
of [3j Lemma 4.2], see also [HI Lemma 4.4]. 

Lemma 5.3. Assume that T = (T(t)) t >o is a strongly continuous contraction 
semigroup on C{&) and that H : C{ff) — > C{ff) is such that for some constants 
a > and N £ IN we have (H(u + v),<p) < a{l + \\v\\ 00 ) 2N+1 for all u, v £ C(0) 
and <p £ <9||u||oo. If for some r > 0, uq £ C{ff) and two continuous functions 
u, v : [0, r) —¥ C(0) we have 

u(t) = T(t)ua+ f T(t-s)F(u(s) + v(s))ds Vt £ [0,r), 
Jo 

Then 

||«(*)||oo< ||«o||oo+ / a(l + |Ks)||oo) 2Ar+1 <fe Vt£[0,r). 
Jo 

Proof of Theorem \5. 6 A Throughout, we fix a positive time T > 0. It suffices to 
prove existence of solutions defined for times t £ [0, T]. Indeed, by the results of 
Section [3] and the Yamada-Watanabe theorem, any solution exist strongly and any 
two solutions agree as long as they are both defined. Thus, solutions defined on 
bounded time intervals can be glued together to get a solution defined for all times. 
We now proceed in two steps. 

Step 1 Construction of a maximal solution for equation (|1.2j) . 
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To begin with, note that by our additional assumptions (1) and (2), the operators 
Ki and Gi are uniformly bounded, i.e. there exists a constant M > such that 

\\Hl( u i> ■ ■ ■ i u r)\\c(&) <M and \\Gi{u)\\^(, 2 ,LP{e)) < M 

for all u\, . . . ,u r G C(ff) resp. u € C(ff). We now also approximate the functions 
Hi with uniformly bounded functions. To that end, we set 

h {n) (x s) := [ ^( x ' s )' if ' s l - n 

1 ' ' 1 hi(x,n ■ sgns) otherwise. 

for n G I and I = 1, . . . ,r. We define H^ : C[0) -> C(^) by [if { (n) (u)] (x) := 
/i, (n) (x,u(x)) for u G C(^) and set F\ n) := K t + H^ n) . Then, for fixed n G IN, the 
map F<") := (i^,. . . , F r (n) ) : C(^) r -> C(^) r is uniformly bounded. 

It follows from [3j Theorem 4.5] that there exists a solution u n = (u± , . . . , Ur) 
of the stochastic equation with coefficients A, F^ n ' and G and initial datum £. By 
pathwise uniqueness these solutions exist strongly. Thus we may assume that the 
processes u n are defined on a common stochastic basis (f2,S,F,P) and that they 
are solutions with respect to a common H- Wiener process Wh ■ Define 

n 

p n :=inf{ie [0,T] : £ |K (n) (*)lloc > n} 

with the convention that inf = T. 

Since F^ n '(u) = F^ n+1 ^(u) for u £ E with \\u\\e < ^, we can repeat the 
arguments in the proof of Theorem 13.11 to prove that u n (t) = u n+ \(t) for all 
< t < p n A This implies in particular that p„ < p n +i almost surely. Thus, 

putting p := sup n/ o„, we may define a maximal solution of our original problem 
dHU) by (w(*))ie[o,p)i wh ere 

:= u„(i) for f G [0,p n ] . 

Step 2 We show that the maximal solution (u(i))te[o,p) is globally defined. 
We prove that for p > 2 we have 

(5.1) sup [|«n||ii>(n;Cf(ro,T];C(^) r ) < °°- 

By [HI Corollary 2.6] this will then imply that p = T and that u is an element of 
L p (rt; C([0, T}; C(ff) r )). It now easily follows that u is a solution of our problem 
ifO)! on the interval [0,T]. 

To prove (|5.1I) . we define the stochastic processes and io[ n ' for n G IN and 
/ = 1, . . .r by 

^ (n) := S t * K^) and W; (n) (i) := S t o l (?,(«,) . 

Since 77; and G; are uniformly bounded it follows from Lemma l5.1l that the processes 
v\ n) and w\ n) are uniformly bounded in LP(fl; C([O,T\;C(0)), say by C p . 

Now put yf 1 ' :— ui — — tuj . Since ti„ is a mild solution of the approximate 
problem, it follows that 

y[ n \t) = S l m + f S t (t s)Hl n \yl n \s)) ds + v[ n \t) + «f>(t) 
Jo 

almost surely for all t G [0,T]. Note that flj satisfies assumption (Fl) with the 
same constants as hi . Consequently, Lemma 12.31 and Lemma 15.31 yield for p > 2 
that 

E||i/ (n) || p _ 

" y l »C([0,T];C(&)) 
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< E sup (|Mloo+ / a / (l + \\v^\s)+w { l n) (s)\\ 00 ) 2N '+ 1 ds 



v 



te[o,T] v Jo 




< 1 + 2C (2Nl+1)p =: M l 



It thus follows that 




C([0,T];C(&)) 



C([0,T];C(^)) 



Summing over / = 1, . . . , r, (|5.ip is proved. 



□ 



Note that assumption (F2) was not used in the proof of Theorem 15.21 That 
assumption is used to get rid of the additional boundedness assumption on the G; 
and Ki in Theorem 15.21 The main tool we use is the following Lemma, which can 
be found in p~4j Lemma 4.8]. 

Lemma 5.4. Assume that T = (T(i)) t > is a strongly continuous contraction 
semigroup on C{&) and that H : C{ff) — > C{&) is such that for some constants 



a,b>0 andN eJN we have (H(u + v) - H(v), ip) < a(l + \\v\\oo) 2N+1 - b\\u\\ 
for all u,ve C{0) and <p G 0||u||oo. If u,v € C([0, T]\C{ff)) satisfy 



Theorem 5.5. Assume (A), (F) and (G). Let (fi, £, F, P) be a stochastic basis 
on which an H -cylindrical Wiener process Wh with respect to F is defined. Then 
for every initial datum £ = (£1, . . . , £ r ) € L (Q, =^o>lP) C(0) r ) there exists a unique 
solution u = (iti, . . . , u r ) of equation (|1.2p . 

Proof. To proof existence of solutions we approximate the coefficients Gi and Ki 
with bounded coefficients. To that end, we define k\ n ^ by setting k\ n \x,s) :— 
ki(x,s) if ||s||i < n and k\ n \x,s) := ki(x, n\\s\\i 1 s) otherwise. The associated 

(n) (n) 

composition operator is denoted by K{ . The functions g\ and the associated 
composition operator are defined analogously. Clearly, and Gj are uniformly 
bounded. Moreover, they are of linear growth with constants independent of n. 

Setting F[ n) := Hi + K\ n) and F< n ) := (F 1 (n) , . . . , F r (n) ), it follows from Theorem 
15.21 that there exists a solution u n — (u^, . . . , it£ ) of equation with coef- 

ficients A,FW and <?( n ) := (G^\ . . . , Gf' l) ). Note that in Theorem UJ we have 
proved existence of solutions only for deterministic initial data. However, as we have 
already remarked in Section 12.51 given pathwise uniqueness we obtain existence of 
solutions for arbitrary J^o-measurable initial data £. 

Proceeding as in the proof of Theorem l5.2[ we can glue together the solutions u n 
to a maximal solution of our original problem and finish the proof of existence by 
proving that the approximative solutions u n for initial data £ € L p (fl, J^q, P; E) are 
uniformly bounded in L p (fi; G([0, T]; E)). This will then prove existence of solutions 
for £ G L P (Q, #q, P; E) . The abstract results of [TJ] will then yield existence of 
solutions for ^-measurable initial data £. 

To prove such a bound, we introduce the processes 




then 




■v. 



Si*K['(u n ) and w\ 



(n) 



:=S,o ; G< n >(u, (n) ). 
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Since the maps K^' and Gv are of linear growth with constants independent of 
n, it follows from Lemma |5. II that for certain constants d\,di and a constant C{T) 
which converges to zero as T — > we have 

e II^ i) Hc ([ o,t ]; cw) * C (W + ^II^ P ,W 

and 

E II^ 3 H P c ([ o,t ]; ^)) ^ C ^ + ^ E E K% M ,c(*) 

Now we put yi := u\ — vj n ^ — wj — S/(-)£/. Since m„ is a mild solution of the 
approximate problem, 

y\ n) (t) = / $ (t - fl )£r, [yl n) (s) + v\ n) (t) + w\ n) (a) + S, (atf,] cfe . 
Jo 

Consequently, by Lemma |5.4[ 

Fll7/ (n) ll p - 

< i + E||t,w+^ +^(^r c([0 , r];C( ^ 

S 1 + E N + ^ Wcao^cm + E H ( *° llc ([ o,T ]; aw) 

Combining with the estimates for v[ n ^ and w\ and summing over I — 1, . . . , r, we 
obtain 

E KIIc([o,t ];B ) < + Ci(T)E||&|[l: + C 2 (T)EK||* ([Q T];B) 
for certain constants Cq,C\(T) and C2(T), where C%(T),C2(T) converge to as 
1 1 0. Thus, if we choose To so small that C^Ib) < 1, we obtain 

1 

l C([0,To];C(«r) - 1 _ C£(T ) ' 

Thus we have obtained the uniform L p -bound for small time intervals. Solving with 
initial datum £ := it(Tb) and with cylindrical Wiener process Wn{t) '■= Wnit+To) — 
Wh(Tq), we obtain solutions u n . By pathwise uniqueness, u n {-) = u n {- + Tq). The 
same arguments as above yields boundedness of u n in L p (il; C([0, To]; E)), hence 
boundedness of u n in L p (f2; C([0, 2Tb]; i?)). Inductively, we obtain boundedness in 
L p (ft; C([0, T\;E)) for all T > 0. □ 

6. An example 

In this section, apply our results to specific stochastic reaction diffusion system 
by presenting a typical example for the reaction term /. More precisely, we consider 
the situation where r = 2, rename the variables u :— u\ and v = V2 and consider 
the function / : R 2 — » R 2 , defined by 

fi(u, v) := u — u 3 + v and f2(u, v) := au — bv 

where a and b are positive constants. This choice for the reaction term goes back 
to the fundamental work of Fitzhugh [5] and Nagumo, Arimolo and Yoshizawa 
[22] . Reaction-diffusion equations involving this reaction term are generic excitable 
systems and appear frequently in applications in chemistry and biology, see e.g. 

This nonlinearity satisfies our assumptions (F). Indeed, setting hi(s) :— s — s 3 
and k\(u,v) := v, then clearly k\ is Lipschitz continuous and of linear growth. 
It follows from Example 12.21 that /i satisfies (Fl) and (F2). Moreover, can put 
/i 2 = 0, which trivially satisfies (Fl) and (F2), and h-i{u,v) := u — v. We should 
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also note that the reaction term / is quasi positive, as /i(0, v) — v > for v > 
and fe(u, 0) = au > for u > 0. 

Thus, our results immediately yield the following 

Theorem 6.1. Let & be a bounded Lipschitz domain in R d , ai,a 2 € L°°(0\ R dxd ) 
be symmetric and uniformly elliptic. Moreover, let R\,R2 be Hilbert- Schmidt op- 
erators on L 2 (0) which are diagonalized by orthonormal bases (e^) resp. (e^) such 
for all k and J2\\ R i e k\\'L,J2\\ R 2ek\\ 2 co < oo. Finally, let 



R be locally ^-Holder continuous and of linear growth. Let (O, S, F, P) 



that ek,e~k 
gx,g 2 : R - 

be a stochastic basis on which two independent L 2 \0) -cylindrical Wiener processes 
W\ and W 2 are defined. 

Then for every £1, £2 G L°(Q, P; C(0)), there exist a pathwise unique solution 
(u, v) of the stochastic reaction diffusion system 



du(t) 
dv{t) 

du(t) 



u(0) 
u(0) 



dvjt) 



div (aiVw(i)) + tt(i) - u(t) 3 
div (fl 2 Vt)(t)) + au(i) - bv(t 
on 9^ 

6 
6- 



+ v(t)]dt + g 1 (u(t))R 1 dW 1 (t) 
)]dt + g 2 (u(t))R2dW2(t) 



If£i,£,2 ore almost surely positive and gi(0) 
almost surely positive for all times. 



52(0) = 0, then the solutions u, v are 
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